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ABSTRACT 

Far-ultraviolet spectra obtained with FUSE show that the strong CI/A1036 in- 
terstellar absorption-line is essentially black in five of the UV-brightest local starburst 
galaxies. Since the opacity of the neutral ISM below the Lyman-edge will be signifi- 
cantly larger than in the CII line, these data provide strong constraints on the escape 
of ionizing radiation from these starbursts. Interpreted as a a uniform absorbing slab, 
the implied optical depth at the Lyman edge is huge (tq > 10 2 ). Alternatively, the 
areal covering factor of opaque material is typically > 94%. Thus, the fraction of ion- 
izing stellar photons that escape the ISM of each galaxy is small: our conservative 
estimates typically yield f esc < 6%. Inclusion of extinction due to dust will further de- 
crease f esc . An analogous analysis of the rest-UV spectrum of the star- forming galaxy 
MS 1512 — CB58 at z =2.7 leads to similar constraints on f esc . These new results 
agree with the constraints provided by direct observations below the Lyman edge in a 
few other local starbursts. However, they differ from the recently reported properties 
of star-forming galaxies at z > 3. We assess the idea that the strong galactic winds 
seen in many powerful starbursts clear channels through their neutral ISM. We show 
empirically that such outflows may be a necessary - but not sufficient - part of the pro- 
cess for creating a relatively porous ISM. We note that observations will soon document 
the cosmic evolution in the contribution of star-forming galaxies to the metagalactic 
ionizing background, with important implications for the evolution of the IGM. 
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mation 

1. Introduction 

The intergalactic medium (IGM) contains the bulk of the baryons in the universe (e.g., 
Fukugita, Hogan, & Peebles 1998). Determining the source and strength of the metagalactic 
ionizing radiation field and documenting its cosmic evolution is crucial to understanding the fun- 
damental properties of the IGM at both low- and high-redshift. The two prime candidates for 
producing the background are QSOs and star-forming galaxies. While the contribution to the ion- 
izing background from QSOs can be estimated with reasonable accuracy, considerably less is known 
about the contribution from galaxies. QSOs alone appear inadequate to produce the inferred back- 
ground (e.g, Madau, Haardt, & Rees 1999), especially at z > 3 where their co-moving space density 
declines steeply with increasing redshift (Fan et al. 2000). Moreover, there are only rather indirect 
constraints on the contribution of galaxies to the ionizing background in the low-redshift universe 
(e.g., Giallongo, Fontana, & Madau 1997; Devriendt et al. 1998; Shull et al. 1999). 

The cosmic history of the formation of massive (ionizing) stars has now been determined 
in broad terms (e.g., Madau, Pozzetti, Sz Dickinson 1998; Steidel et al. 1999; Barger, Cowie, 
& Richards 2000). On-average several thousand ionizing photons are produced for each nucleon 
processed into stars, so the galactic contribution to the ionizing background is potentially both 
substantial and calculable. However, the greatest single uncertainty in this calculation is the value 
for f esc - the fraction of the ionizing photons that escape from star-forming galaxies. 

Uncertainties in the value of f esc also map directly into uncertainties in the use of standard 
emission-line diagnostics to infer the basic properties of star-forming galaxies. Most crucially, the 
luminosity of the Hydrogen recombination-lines is a widely utilized measure of the formation rate 
of massive stars in both the local and high-redshift universe (e.g. Gallego et al. 1995; Yan et al. 
1999; Pettini et al. 2001). This approach generally assumes that f esc << 1. However, the smaller- 
than-predicted equivalent widths of the Balmer emission lines in starbursts could imply that this 
assumption is suspect (e.g. Moy et al. 2000; Stasinska et al. 2001). 

An HI column of only 1.6 x 10 17 cm~ 2 produces r = 1 at the Lyman edge. Mean galactic 
gas columns are much larger of course, ranging from ~ 10 21 cm~ 2 in normal galactic disks to 10 24 
cm~ 2 in nuclear starbursts (e.g., Kennicutt 1998). Starburst-like mean gas columns are inevitable 
in the Lyman Break Galaxies (Heckman 2000), given their high star- formation rates per unit area 
(Meurer et al. 1997). The leakage of ionizing radiation out of galaxies must then be determined 
by the topology of the ISM. As such, theory gives us scant guidance, and so direct measurements 
of f esc are required. 

Leitherer et al. (1995) reported the first direct measurements of f esc using the Hopkins Ul- 
traviolet Telescope to observe far-ultraviolet (far-UV) light below the rest-frame Lyman edge in 
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a sample of four local starbursts, and these data were later reanalysed by Hurwitz, Jelinsky, & 
Dixon (1997). The resulting upper limits on f esc were typically 3% to 10%. In a very surprising 
development, Steidel, Pettini, & Adelberger (2001 - hereafter SPA) have reported the detection 
of escaping Lyman continuum in the combined spectrum of 29 Lyman Break Galaxies at a mean 
redshift < z > = 3.4. They estimate that the ratio of f esc at 900 A and 1500 A ranges from 0.5 to 1 
in this composite spectrum (i.e., there is no definite detection of photoelectric opacity due to HI). 
Haehnelt et al. (2001) reach similar conclusions, but point out that star-forming galaxies at only 
slightly smaller redshifts can not be this porous without violating constraints set by the observed 
Hell/HI opacity ratio in the Lya forest. 

It is clearly important to obtain more determinations of f esc , and to use these measurements 
to better understand the physical processes that determine f esc . This latter goal can be best 
accomplished in the local universe, where detailed multiwaveband observations of star-forming 
galaxies can be made. Unfortunately, direct observations of the emerging Lyman continuum require 
observing galaxies at redshifts greater than a few percent, so that the foreground HI opacity of the 
Milky Way is not significant at the relevant wavelengths (Leitherer et al. 1995; Hurwitz, Jelinsky, & 
Dixon 1997). Given the modest sensitivity of far-UV telescopes, this has limited such investigations 
to rather small (and possibly unrepresentative) galaxy samples. 

Against this backdrop, we describe our analysis of new far-UV data obtained with FU SE for 
the UV-brightest local starbursts (section 2). We point out that the strongest interstellar lines that 
trace the HI phase in these starbursts are black (or very nearly black) at line-center. We show that 
this implies low values for f esc (section 3). This not only more than doubles the sample of local 
starbursts with good upper limits on f esc , it illustrates a technique that can potentially be used on 
a much larger sample of galaxies. Indeed, similar arguments can be applied to high-redshift galaxies 
(section 4), and show that f esc is also low in the bright star-forming galaxy MS 1512 — CB58 at z 
=2.7 (Pettini et al. 2000). We discuss these results and their implications in section 5. 

2. Observations 
2.1. The Sample 

We are carrying out a program with the Far Ultraviolet Spectroscopic Explorer (FUSE - 
Moos et al. 2000) to obtain spectra of the six starburst galaxies from the Kinney et al. (1993) 
'Ultraviolet Atlas of Star-Forming Galaxies' having the largest UV flux at 1500 A through the 
IUE 10" x 20" aperture (F\ > 10~ 13 erg cm~ 2 s~ - 1 A" 1 ). These targets span broad ranges in: 1) 
metallicity, from 1/8 (NGC 1705) to 2.5 (M 83) times solar, 2) starburst bolometric luminosity, 
from ~3xl0 8 (NGC 1705) to -2xl0 10 (NGC 3310) L , 3) internal dust-reddening, from Ay ~ 0.0 
(NGC 1705) to 0.6 (NGC 3310 and M 83), and 4) host galaxy properties, from dwarfs (NGC 1705 
and NGC 5253), to irregulars (NGC 4214 and NGC 4449), to spirals (NGC 3310 and M 83). The 
sample therefore spans most of the multidimensional parameter space of starbursts in the local 
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universe (Heckman et al. 1998). The selection on the basis of rest-frame UV brightness makes this 
an appropriate sample to compare to the UV-selected Lyman Break Galaxies at high-redshift. 

2.2. Observational Details 

The observations undertaken to date are summarized in Table 1. Data have been obtained 
for 5 of the 6 targets (only NGC 4449 has not yet been observed). Analyses of the dynamics of 
the ISM in individual galaxies based on these data are being reported elsewhere: Heckman et al. 
(2001a,b), Martin et al. (2001). The large (LWRS, 30" x 30") aperture on FUSE was used for 
the observations of 4 of the 5 targets, while NGC 5253 was observed with the medium (MDRS, 
4" x 20") aperture. The corresponding physical sizes of the projected aperture are kpc-scale in all 
cases but NGC 5253 (Table 1). 

The starburst was centered in the aperture of the LiFl (guiding) channel for each observation 
by the standard guide-star acquisition procedure. For four targets, flux was recorded through 
the LWRS apertures in both long wavelength (LiF, ~ 1000 — 1187 A) channels and both short 
wavelength (SiC, ~ 900— 1100 A) channels. We obtained no useful data in the SiC channels during 
our MDRS observations of NGC5253 (A0460404) due to a loss of mirror alignment. We note 
that while the data do extend slightly below the wavelength of the starburst's Lyman edge in the 
SiCl channel, these data do not usefully constrain f esc . This is because the opacity of foreground 
Galactic HI is still significant at these wavelengths, largely due to the confluence of the high-order 
Lyman series absorption lines (cf. Hurwitz et al. 1997). 

Most of the data for this program were obtained at night, which minimizes O I and N I 
terrestrial airglow contamination in the spectra. The percentage of night to total exposure time 
averaged ~ 65%, except for NGC 1705 for which the ratio was ~ 92%. 

The velocity resolution of the data depends on how the far-UV light illuminates the FUSE 
apertures. The resolution for a filled MDRS aperture (NGC 5253) or for a point-like source in 
the LWRS aperture (NGC 1705) is ~ 30 km s _1 . We estimate a resolution of <50 to <70 km s" 1 
in NGC 3310, M 83, and NGC 4214 based on the widths of the narrowest Galactic (foreground) 
absorption-lines. 

2.3. Data Processing 

The raw time-tagged photon event lists for each exposure were processed with the standard 
FUSE calibration software (C ALFU S Evl.8.5) available at the Johns Hopkins University as of 
January 2001. The lists were screened for valid data, and corrections for geometric distortions, 
spectral motions, and Doppler shifts were applied (see Sahnow et al. 2000). The individual cali- 
brated extracted spectra for each channel were cross-correlated, shifted to remove residual velocity 
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offsets due to image motion in the apertures, and combined to produce a composite spectrum. 
These composite spectra in the 1000-1070 A region were then compared, and any remaining ve- 
locity offsets were removed by referencing the Galactic absorption lines to the appropriate Local 
Standard of Rest velocities. 

Since we are interested in measuring the fluxes in the cores of saturated lines, accurate de- 
terminations of the residual fluxes are essential. We applied the standard detector background 
corrections available within CALFUSE for each observation. These corrections account only for 
the particle event backgrounds (~ 0.75 cnt cm -2 s" 1 ) and do not correct for scattered Lya emission 
or other stray light. These additional sources of background light are negligble in the LiFl data at 
the wavelengths of ~ 1040 A (typically < 5 x 10~ 15 erg cm -2 s^ 1 A -1 , or less than < 2% of the con- 
tinuum flux). A description of the FUSE backgrounds can be found in the FUSE Observers Guide 
(Blair & Andersson 2001) 3 . 

In the following discussion we consider data from the LiFl channel, which has the highest S/N 
at the wavelengths of interest for this study (A ~ 1020 — 1045 A). The data for the other channels 
are consistent with our findings presented below. 

3. Results from the FUSE Data 

In Figure 1 we show the spectral region around the Ly(3, C II A1036.337 and O I A1039.230 
absorption-lines in the five galaxies in our sample. It is immediately apparent that the cores of the 
C II A1036 and Lyman lines are quite black. The residual normalized intensities in the C II line are 
typically < 6% (Table 2). These C II residual intensity estimates account for overlying Galactic 
absorption due to the Milky Way O VI A1037.617 or H 2 R(l) A1037.149 and P(l) A1038.157 lines. 
There are no H2 lines within the starbursts themselves that contribute to the depth of the C II line. 
We modeled the Galactic absorption features using the shorter wavelength member of the O VI 
doublet (A1031.926) or other H 2 R(l) and P(l) lines (typically P(l) AA1003.294, 1014.326 and R(l) 
A1013.435) to get both the velocity extents and expected depths of the overlying lines. We note that 
the calculation of the residual intensity is also corrected for the effect of the maximum plausible 
stellar contribution to the C II line. Model starburst spectra (Gonzalez Delgago, Leitherer, & 
Heckman 1997) show that the stellar photospheric line (which is entirely due to B stars) will depress 
the local continuum by a factor of ~ 1.2 to 2 (primarily depending on the age of the starburst). In 
our sample, the stellar contribution to C II should be largest in the B-star-dominated spectrum of 
NGC 1705. Even in this case the C II line is primarily interstellar, since it (like the other interstellar 
lines) is blueshifted by roughly 40 km s _1 relative to the galaxy systemic velocity (Heckman et al. 
2001, and see also Heckman & Leitherer 1997). 

The O I A1039 line shows behavior ranging from a normalized residual intensity of about 15% 



3 This guide (v3.0) can be found on-line at http://fuse.pha.jhu.edu/support/guide/guide.html. 
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Fig. 1. — FU SE spectra (LiFIA detector segment) of the spectral region containing the interstellar 
Ly(3, C II A1036, and O I A1039 absorption-lines in five UV-selected local starbursts. The lines due 
to the foreground (Galactic) and starburst ISM are indicated, as is the fitted continuum (dotted 
line) to which the residual intensity at the line core is referenced. Note that the Lyj3 and C II A1036 
lines are effectively black in the line core, while the less saturated O I A1039 line shows a range in 
residual intensity. Tick marks indicate the positions (or expected positions) of J = — 3 H2 lines 
in the Milky Way (upper, lighter ticks) and the starbursts (lower, darker ticks). All of the spectra 
shown are a combination of data taken during both orbital day and night. 
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in NGC 1705 to < 6% in NGC 5253. Overlying Galactic absorption is not an issue for the O I 
at the resolution of the FUSE data. Both O I and C II arise in H I regions and contain direct 
information about the H I column density and areal covering factor of the neutral gas. 

We begin by writing the ratio of the optical depth at the Lyman edge to the optical depth at 
line center. Following Spitzer (1977), this is 



where (Nn/Nj) is the ratio of the H I and ground-state ionic columns, b is the normal Doppler 
parameter (cm s' 1 ), f is the oscillator strength, and the wavelength A is measured in cm. Assuming 
that O I and C II are the dominant ionic species of their respective element in the H I phase, this 
can be written as 



where Z gas is the gas-phase C or O abundance in solar units and R = 8.6 for C II A1036 and R = 
50 for O I A1039 (Morton 1991). The important point is that the optical depth at the Lyman edge 
is much larger than the optical depth at line-center for either line for typical values of Z gas and b. 
Corrections for grain-depletion will only increase the implied ratio of tl^/to- Ionization corrections 
are unimportant for O I since it is closely coupled to H I through charge exchange reactions, and 
will only increase the implied ratio of tl v /tq for C II. 

As an example, consider the case of NGC 1705 (Heckman et al. 2001). Based on the neb- 
ular emission- lines, the gas-phase Oxygen abundance is 1/8 solar, and this is consistent with the 
abundances derived in the H I gas for Si, Ar, and Fe from the FUSE spectra. The measured b 
value is 40±10 km s" 1 . The residual intensity at the core of the O/A1039 line implies tq = 1.8 (for 
unit covering factor - see below). Equation 2 above then implies that tl v = 280. This is roughly 
consistent with the HI column we derived from fitting the Lyman series lines: Njji = 1.5x 10 20 
cm -2 . Thus, in the context of a homogeneous absorbing slab, the escaping fraction of ionizing 
photons will be almost identically zero in NGC 1705. In principle, even more severe constraints on 
TL y can be placed using weaker absorption-lines (transitions having smaller values of (Nj/Nu)fX)- 
We do not take this approach, since even the strong O I A1039 line already implies f esc ~0. 

Less stringent (more realistic) constraints on f esc can be placed by adopting a "picket fence" 
model in which the areal covering factor of optically-thick HI is not unity (e.g., ionizing radiation 
can escape through "holes" in the HI). In this case, the upper limit on the residual intensity at 
the center of the C//A1036 line implies that the covering factor must be > 94%, so that the upper 
bound on the escape fraction for ionizing radiation is f esc < 6%. 

We have undertaken this analysis for the other four galaxies in our sample, and report the 
results in Table 2. In all cases, the results are similar to those in NGC 1705, with implied values for 



tlv/to = 4 x lO-^AVA^XV/A) 



(1) 



T Ly /To = RZ-i(b/100 km/s) 



(2) 
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f esc that are almost identically zero (for the homogeneous slab model) and <6% for the picket-fence 
model. The values listed in Table 2 should be considered stringent upper limits for the following 
reasons: 1) We list only the limits derived from the picket fence model. 2) We measured f esc 
at velocities where the residual flux was greatest after accounting for possible contamination by 
other absorption features. 3) Scattered light may contribute at low levels (< 2%) and only serves 
to strengthen the limit on f esc . 4) The FUSE line spread function (LSF) likely consists of two 
components: a narrow component accounting for ~ 80% of the LSF area and a broad component 
accounting for ~ 20% and having ~ 2 — 3 times the width of the narrow component. While the 
exact details of the LSF are presently unknown, such a model more accurately describes the shapes 
of the cores of the Galactic H I lines toward nearby stars than a single component LSF. As a result, 
it is possible that some of the observed flux in the line cores is simply redistributed light from 
nearby continuum regions within ~ 1 A of the line. 5) The continuum estimates for the interstellar 
lines shown in Figure 1 are conservative estimates and are unlikely to be lower than those shown. 
6) The upper limits on f esc can be pushed to even smaller values by including the opacity due to 
dust (that is, our technique measures only the contribution of the photoelectric opacity of H I). 
Application of the empirical starburst dust attenuation law (Calzetti 1997; Meurer, Heckman, &; 
Calzetti 1999) to the IUE spectra of these galaxies implies an extinction due to dust at 1500A 
ranging from negligible (NGC 1705) to ~ 3 magnitudes in M 83 and NGC 3310. The mean dust 
opacity in the Lyman continuum (1 to 4 Ryd) should be larger still (Mathis 1990). 

We conclude that our UV-selected starbursts are quite opaque to their own ionizing radiation. 
This confirms the results of Leitherer et al. (1995) and Hurwitz, Jelinsky, & Dixon (1997), but 
generalizes and strengthens them. 

4. Galaxies at High-Redshift 

The technique described above can be applied to galaxies at high-redshift. As an application, 
we consider the star-forming galaxy MS 1512 — CB58 at z =2.7 (Yee et al. 1996). Thanks to 
strong gravitational lensing, its flux has been amplified by a factor of ~ 30 (Seitz et al. 1998). 
Thus, the rest-frame UV spectrum discussed by Pettini et al. (2000) is by far the best of its kind. 
We display the spectral region from Lya to C//A1335 in Figure 2. The quoted spectral resolution 
in these data is ~200 km s~ x FWHM, while the interstellar lines are much broader (500 to 700 km 
s _1 ) and very deep. It is instructive to note that the value for /A (equation 1) for the CT/A1335 
(O/A1302) line is 0.1 (0.8) dex larger than for the C 11X1036 (O7A1039) line in our FUSE spectra. 
Thus, the spectrum of MS 1512 — CB58 can be straightforwardly compared to the local starbursts. 

The measured residual intensites in the line cores in the MS 1512 — CB58 spectrum are: Lya 
(1%), Sil 1X1260 (7%), O/A1302 (6%), and C 11X1335 (2%). To apply equations 1) and 2) above, 
we adopt a gas-phase metallicity of 1/4 solar (Pettini et al. 2000; Leitherer et al. 2001; Tremonti 
et al. 2001) and take b = 300 km s _1 based on the line-widths reported by Pettini et al. (2000). 
The implied optical depth at the Lyman edge assuming a uniform slab is tl v ~ 300. Pettini et al. 
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Fig. 2. — A portion of the Keck LRIS spectrum of the gravitationally-lensed star-forming galaxy 
MS 1512 - CB58 at z =2.7 published by Pettini et al. (2000). The interstellar Lya, Si II A1260, 
O I A1302, Si II A1304, and C II A1335 absortion-lines are indicated. The stellar continuum in the 
spectrum has been normalized to unit intensity by fitting it with a low-order polynomial. Note 
that the interstellar lines are nearly black, and much broader than the instrumental resolution of 
0.8 A (rest-frame). 
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(2000) derive an HI column of 7.5 xlO 20 cmT 2 from the damped Lya line, yielding an even larger 
value for tl v for a uniform slab. Alternatively, the residual intensity at the core of the C//A1335 
line implies an areal covering factor of optically-thick gas of 98%. The implied value for f esc ranges 
from ~0 to 2%. These values are consistent with the results for the local starbursts in section 3 
above. 

The situation is quite different in the faint Lyman Break Galaxies at z > 3. As SPA point 
out, the interstellar lines in their composite spectrum are far from black (the apparent residual 
intensities are roughly 70% at line center). The composite spectrum of Lyman Break Galaxies at 
< z > = 3.0 in the Hubble Deep Field is similar (Lowenthal et al. 1997). It is possible in these 
cases that spectra with higher spectral resolution and signal-to-noise would show that the lines are 
in fact nearly black over narrow ranges in velocity. However, the existing data are consistent with 
the conclusion of SPA that the escape fraction of ionizing radiation is large in these galaxies. 

5. Discussion 

The results we have reported for five of the six UV-brightest local starbursts (NGC 1705, 
NGC3310, NGC 4214, NGC 5253, and M 83) implies that UV-bright starburst galaxies today 
are highly opaque to their own ionizing radiation. This agrees with results for a small sample of 
more distant starbursts described by Leitherer et al. (1995) and Hurwitz et al. (1997). These 
results in turn are consistent with the estimate by Shull et al. (1999) that the strength of the local 
metagalactic ionizing background implies that the fraction of ionizing photons that escape from 
star- forming galaxies can not exceed a global average value of about 5%. 

These results are in striking contrast to the detection of significant escaping Lyman continuum 
radiation from Lyman Break Galaxies at z ~ 3 to 4 by SPA. Our analysis shows that the star- 
forming galaxy M51512 — C-B58 at z =2.7 is similar to the local starbursts: it too must be highly 
opaque below the Lyman break. Taking the current results at face value, we are left with a puzzle. 
Why is the topology/structure of the neutral interstellar medium apparently so different between 
the galaxies studied by SPA and those in the present paper? 

In some ways this is a surprising result, because local starbursts and Lyman Break Galaxies 
are remarkably similar in most respects (e.g., Heckman 2000). Meurer et al. (1997) showed that 
both have the same bolometric surface-brightnesses, and therefore the same star formation rates 
per unit area. This implies that the basic physical properties of the ISM should be similar. The 
surface-mass-densities in the stars (£*) and the interstellar gas (S 9 ) that fuels the star-formation 
must be similar. The radiant energy density in the ISM must be similar (~10 3 times higher than 
in the Milky Way) , as must the rate of mechanical energy deposition (supernova heating) per unit 
area or volume. Simple considerations of hydrostatic equilibrium imply correspondingly high total 
pressures in the ISM in both: P ~ GX^E^, or P/k ~ few x 10 7 K cm -3 . This is several thousand 
times the value in the local ISM in the Milky Way. Finally, the characteristic dynamical time in 



-li- 



the ISM will be short in both cases: t dyn ~ (Gp) _1/2 ~ (GY ltot H)~ 1 / 2 ~ few Myr (where # ~ 10 2 
pc is the thickness of the disk). Perhaps the most striking difference is the size-scale: starbursts 
in the local universe are generally circum-nuclear events (sizes of 10 2 to 10 3 pc) imbedded in a 
much larger galaxy, while the intense star-formation in the Lyman Break Galaxies appears to be a 
galaxy- wide event (sizes of 10 3 to 10 4 pc). This might lead to fundamental differences in the global 
properties of the ISM. 

Such generic arguments aside, is there any empirical evidence that the structure or dynamics 
in the ISM in our specific sample of local starbursts differs from the Lyman Break Galaxies? We 
consider two possibilities. 

First, our FUSE sample galaxies have considerably lower luminosities (star-formation rates) 
than the typical Lyman Break Galaxies: logL^i ~ 8.5 to 10.3 L Q vs. 10.7 to 11.7 L© (Adelberger & 
Steidel 2000). It seems physically plausible that the ISM might become more porous in starbursts 
with higher star- formation (and energy deposition) rates (cf. Lehnert & Heckman 1996). We have 
tested this idea by considering data on f esc for more powerful local starbursts. Two of the four 
starbursts observed below the Lyman break by Leitherer et al. (1995) are powerful systems with 
Lbol ~ 10 11 L Q . Both are opaque in the Lyman continuum: f esc < 2% and < 4% for Mrk 496 and 
IRAS 08339+6517 respectively (Hurwitz et al. 1997). We have also searched the HST archives for 
UV spectra of powerful local starbursts taken with a spectral resolution comparable to our FUSE 
data and with a projected aperture size (~kpc-scale) comparable to our FUSE data. There are 
three such starbursts in the sample investigated by Kunth et al. (1998), who used the GHRS echelle 
mode to observe the Lya, O/A1302, and Szi7A1304 lines with a resolution of ~ 20 km s -1 . These 
starbursts have luminosities (L ) of logL hol = 10.4 (ESO 400-G043), 11.0 (IRAS 08339+6517), 
and 11.2 (ESO 350-IG038). These interstellar lines are black in two cases (ESO 400-G043 and 
ESO 350-IG038), and have a residual intensity of 60 to 70% in the third case (IRAS 08339+6517). 
While a large value for f esc is therefore possible in this galaxy, direct observations below the Lyman 
break actually show that f esc < 4% in this case (Hurwitz et al. 1997). 

A second possibility is that large values for f esc are created by supernova-driven galactic 
outflows (superwinds) that clear out channels through which ionizing radiation can escape (SPA). 
The tracers of superwinds in Lyman Break Galaxies are their blueshifted interstellar absorption- 
lines (Franx et al. 1997; Pettini et al. 1998; 2000; 2001). Indeed, the composite spectrum of the 29 
Lyman Break Galaxies analyzed by SPA shows the clear signature of these outflows: the interstellar 
absorption-lines are blueshifted by 400 to 500 km s _1 with respect to the Lya emission-line. This 
outflow signature is also seen in the composite spectrum formed from the sum of the spectra of 12 
Lyman Break Galaxies in the HDF at z ~ 3 (Lowenthal et al. 1997; Franx et al. 1997). As noted 
above, the relative shallowness of the interstellar lines in both these composite spectra is consistent 
with (but does not demand) a large value for f esc . 

Similar outflows are seen in the UV spectra of local starbursts (Kunth et al. 1998; Gonzalez 
Delgado et al. 1998), and Heckman et al. (2000) have shown that high- velocity outflows of neutral 
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gas are common in powerful (L^ > 10 11 L Q ) local starbursts. In contrast, our FUSE data reveal 
high- velocity (Av > 10 2 km s™ 1 ) outflows in the neutral phase of the ISM in only one of our five 
galaxies: NGC 3310, in which the mean outflow velocity is about -200 km s _1 (see Figure 1). In the 
other galaxies the neutral-phase absorption-lines arise in relatively quiescent gas near the galaxy 
systemic velocity (Figure 1; Martin et al. 2001; Heckman et al. 2001). Based on the ISM dynamics, 
it is perhaps not surprising that NGC 1705, NGC 4214, NGC 5253, and M 83 are opaque to their 
Lyman continuum radiation. 

On the other hand, Gonzalez Delgado et al. (1998) have shown that outflows at several hundred 
km s _1 are present in the neutral gas in at least three of the four starbursts observed below the 
Lyman break by Leitherer et al. (1995). In all cases the galaxy is opaque below the Lyman break, 
with f esc typically < 3 to 10%. Of the three powerful starbursts observed by Kunth et al. (1998), 
two (ESO 400-G043 and ESO 350-IG038) show black interstellar lines and outflows at (v — v sys ) 
= -225 and -58 km s _1 respectively. In the third (IRAS 08339+6517), the lines have a residual 
intensity of 60 to 70%, and are strongly blueshifted by ~-500 km s" 1 (Gonzalez Delgado et al. 
1998). 

We conclude that the currently avalilable data do not demonstate that galactic winds inevitably 
produce large values of f esc in local starbursts. Such outflows appear to be a necessary - but not 
sufficient - part of the process that creates an ISM that is porous to ionizing radiation. 

In the near future we will use FUSE to observe a sample of more powerful starbursts (L^i 
= 4xl0 10 to 4xl0 n Lq). This will allow us to probe the Lyman continuum opacity in more 
starbursts with luminosities (and star formation rates) that are similar to those of typical Lyman 
Break Galaxies. At the same time, the extension of the investigation at high-redshift to a wider 
range of Lyman Break Galaxies can help clarify the situation. In particular, it will be important 
to test the high values for f esc by obtaining data on both the Balmer recombination-lines and the 
Lyman continuum in the same set of galaxies (to provide a consistency check) . Finally, the Galex 
mission (Martin et al. 1999) will directly measure the escaping flux below the Lyman break for the 
population of star-forming galaxies in the redshift range z ~ 0.4 to 2 (the epoch that apparently 
dominates the cosmic history of star- formation). Thus, over the next few years, we should be 
able to make the first good direct measurement of the contribution of star-forming galaxies to 
the metagalactic ionizing background as a function of cosmic epoch. This will have important 
implications for the history of the IGM. 

We thank the members of the FUSE team for providing this superb facility to the astronomical 
community. We thank Max Pettini for enlightening discussions, and for kindly providing us a copy 
of the composite Lyman Break Galaxy spectrum published by SPA. This work was supported in 
part by NASA grants NAG5-6400 and NAG5-9012. KRS acknowledges support from NASA Long 
Term Space Astrophysics grant NAG5-3485. 
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Table 1. FUSE Starburst Observations 



Galaxy 


I 

(°) 


b 

(°) 


F : 

1050 A 
(erg cm -2 s _1 


A- 1 ) 


Dataset 


Obs. Date 
(y/m/d) 


rp a 

± exp 

(ksec) 


Aper. b 


L c 
(kpc) 


NGC 1705 


261.08 


-38.74 


~ 4.2 x 10" 


-13 


A0460102 


2000/02/04 


7.7 


LWRS 


0.9x0.9 












A0460103 


2000/02/05 


13.6 


LWRS 




NGC 3310 


156.01 


+54.06 


~ 2.6 x 10" 


-13 


A0460201 


2000/05/05 


27.1 


LWRS 


2.56x2.56 


NGC 4214 


160.25 


+78.08 


~ 1.6 x 10" 


-13 


A0460303 


2000/05/12 


20.7 


LWRS 


0.50x0.50 


M 83 


314.58 


+31.97 


~ 3.0 x 10" 


-13 


A0460505 


2000/07/06 


26.5 


LWRS 


0.55x0.55 


NGC 5253 


314.86 


+30.11 


~ 8.9 x 10" 


-14 


A0460404 


2000/07/07 


27.4 


MDRS 


0.40x0.08 d 



a On-target exposure time in the LiFl channel. 

b Aperture sizes. MDRS: 4x20 arcsec; LWRS: 30x30 arcsec. 

c Physical size of projected FUSE aperture. 

d Aperture position angle = 195.87°. No useful SiC data obtained for this MDRS observation. 
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Table 2. Residual Intensities in the Core of the C II 1036.337 A Line 



Galaxy 


V.sys 

(km s _1 ) 


fesc(C II) 
(%) 


Notes on f esc and nearby absorption 


NGC 1705 


569 


< 5.8 


Weak Galactic high velocity VI A1037.617 in red wing of profile. 
No H2 absorption at nearby velocities. 


NGC 3310 


997 


< 6 - 30 


Value of f esc increases along redward wing of absorption. 


NGC 4214 


298 


< 4.0 


Weak Galactic H 2 R(l) A1037.149 on blue side of profile. 


M 83 


517 


< 6.0 


Weak Galactic H 2 P(l) A1038.157 on red side of profile. 
Galactic O VI A1031.926 on blue side of profile. 


NGC 5253 


405 


< 6.0 


Weak Galactic H 2 P(l) A1038.157 on red s ide of profile. 
Galactic O VI A1031.926 on blue side of profile. 



